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Abstract 

We consider the ordinary differential equation 

x 2 u" = axu' + bu — c(u' — l) 2 , x € (0, xq), 

with o£K,()(!l, c>0 and the singular initial condition it(0) = 0, which in financial 
economics describes optimal disposal of an asset in a market with liquidity effects. It 
is shown in the paper that if a + b < then no solutions exist, whereas if a + b > then 
there are infinitely many solutions with indistinguishable asymptotics near 0. More- 
over, it is proved that in the latter case there is precisely one solution u corresponding 
to the choice xq = oo which is such that < u(x) < x for all x > 0, and that this 
solution is strictly increasing and concave. 
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Introduction 

The paper is concerned with solutions of the problem 



* Corresponding author: brunovsky@fmph.uniba.sk; tel:+421 265425741; fax:+421 265412305 

* ales, cerny. 1 @city. ac.uk 

* michael.winkler@uni-due.de 



84248 Bratislava, Slovakia, 

Ales Cernyf 
Cass Business School, City University London, 
106 Bunhill Row, London EC1Y 8TZ, UK 

Michael Winkler* 
Institut fur Mathematik, Universitat Paderborn, 
33098 Paderborn, Germany 



x 2 u" = axu' + bu — c(u' — l) 2 , 
u(0) = 0, 
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using the scaling 



where a and b are real numbers and c, x$ > 0. By a solution of (jO.ip we mean a function 
u £ C°([0,xo]) n C 2 ((0,xo)) which satisfies (jO.ip in the classical pointwise sense. 

Equation (jO.ip arises in the study of a specific stochastic optimization similar to the 
classical LQ problem. The equation is singular at x = which in itself is not particularly 
noteworthy, since stochastic LQ problems with geometric Brownian state variable invari- 
ably give rise to nonlinear singular ODEs/PDEs of the type seen in (jO.ip and in (|0.3p below, 
see for example [7]. Our problem derives its rich structure from the fact that the initial 
condition (|0.2p . too, refers to the singular point x = 0. This, as we demonstrate below, 
poses certain technical obstacles in establishing existence and, more importantly, gives rise 
to infinitely many solutions with indistinguishable asymptotics near zero (Corollary 13. 4p . 

As was already highlighted, the ODE (10. ip . (|0.2p is not artificial, rather it stems from a 
well-defined optimization problem in financial economics. Specifically, it is obtained from 
the PDE 

f \y 2 a 2 w yy + Xywy + r*zw z - pw + fa ~ 4 ^ } = 0, , Q ^ 

\ io(y,0)=0, 

y 2 z 
w(y, z) = — u(x), x = r]—. 

V V 

The PDE (|0.3p in turn represents the Hamilton-Jacobi-Bellman equation for the opti- 
mal value function w of the following dynamic optimization problem: 

/ r T(z=0) \ 

w(y(0), z(0)) = max E I J e~ ps f(y(s), z(s))(y(s) - V f(y(s), z(s)))ds 1 , (0.4) 
subject to 

dy(t) = \y(t)dt + ay(t)dB(t), (0.5) 

dz(t) = (r*z(t)-f(y(t),z(t)))dt, (0.6) 

over the controls / : M + x M + i— >■ R, T(z = 0) being the first arrival time at z = and 
B(t) the standard Wiener process. 

Optimization (|0.4p with dynamics given by (|0.5l 10. 6p models optimal liquidation of a 
large quantity of an asset whose market price is adversely affected by its ongoing sale. In 
this context z(0) represents the quantity of the asset yet to be sold, y(0) is the prevailing 
price and w captures the expected revenue of an optimal sale of quantity z conditional on 
the current price being y. For more details we refer the reader to [2]. 

The problem of existence is the first principal subject of this paper. It is shown that 
for a + b > the problem (jO.ip , (|0.2p has a continuum of local solutions and at least one 
global solution bounded between and x (Section 1), while for a + b < no solutions exist 
(Section 2). 

If one admits the possibility that there are multiple solutions to (jO.ip , (|0.2p , one imme- 
diately has to deal with the additional challenge of identifying "the" right solution relevant 
to the associated optimization problem. The economic nature of the optimization (|0.4t - 
10. 6p strongly suggests that the relevant solution of (jO.ip . (10. 2p should be increasing (larger 
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amount of asset means larger revenue) but concave (decreasing returns to scale, since 
larger volume of sales has greater adverse effect on the sale price of the asset). However, 
there is no indication in the form of equation (jO.ip that a solution with these properties 
should exist in the first place. In Section 3 we thus analyze monotonicity and convexity 
properties of a solution bounded between and x (Proposition I3.5p . the upper bound 
corresponding to an immediate sale of the entire stock of the asset without any adverse 
price effect. 

In Section 4 we address the question of global uniqueness. We show that there is 
exactly one solution on M + which remains bounded between and x, and this solution 
is necessarily increasing and concave (Proposition 14. Finally, in Section 5 we examine 
finer aspects of local non-uniqueness. In the current paper we focus on the intricacies of 
the initial value problem (jO.ip . (|0.2p . The implication of the results for the underlying 
optimal control problem is a delicate issue left to further research. 

A paper similar in spirit to ours is [3] . It studies a specific second order equation with 
a singularity at 0, arising in the theory of general relativity. As in our case, the existence 
of infinitely many solutions is established by the method of upper and lower solutions and 
then the properties of the set of solutions are studied. 

As far as the local existence is concerned, Liang in [3] carried out a systematic study 
of second order singular initial value problems of the form 



where F is a continuous function and the initial conditions satisfy F(0,u(0),u' (0)) = 0. 
The key quantity in this study is 7 := -^tF(0, u(0), it'(0)). It is shown that for 7 < local 
uniqueness holds, while for 7 > solutions become unique only after the asymptotics of 
v! have been fixed to the order 2 7 near x = 0. The case 7 = is not treated. Each 
solution has an asymptotic expansion in powers of x and x 7 (provided 7 is not an integer) , 
and asymptotic expansion of is obtained by differentiating n-times the asymptotic 
expansion for u. 

In contrast, we study a specific singular IVP from a wider class 



with a = 2, u(0) = 0, u'(0) = 1 and F(x,u,u') := axu' + bu - c(v! - l) 2 . Like in [I] 
our ODE arises from a self-similar solution of a PDE. However, we deal with a borderline 
case where ^F(p,u(0),u' (0)) = 0. As a result, standard blow-up techniques are not 
productive and we have to resort to the method of sub-supersolutions. 

Finally, we remark that it is not uncommon for HJB equations associated with stochas- 
tic optimization to exhibit multiple solutions. The meaningful solution then has to be 
selected by employing additional criteria. In the case of LQ problems the relevant solu- 
tion is identified as the maximal/minimal one. In other cases the optimal solution can be 
singled out as the unique viscosity solution of the HJB equation, cf. p] . In our case these 
criteria do not seem to be helpful. Rather, the significant solution is uniquely determined 
by its global monotonicity and concavity properties. 




(0.7) 



x 



—F(x,u,u') 
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1 Existence for a + b > 

In what follows existence will be shown in a rather direct way by constructing an ordered 
pair of a sub- and a supersolutions to (jO.ip . and applying a standard existence result for 
elliptic boundary- value problems which will provide a solution lying in between, cf. [5], 
[6]. Due to the singularity in the ODE (jO.ip . however, an approximation procedure will 
be involved in the proof. We isolate technical arguments in the following proposition. 

Proposition 1.1. Let < x\ < %2 and suppose that there exist u,u G C°([xi,X2]) Pi 
C 2 ((xi, X2)) such that 

u < u in (x±, X2); (1-1) 
£u < in (x±, X2); (1.2) 
£u > in (xi,X2); (1.3) 

the operator £ being defined according to 

£u := -x 2 u" + axv! + bu- c{v! - l) 2 (1.4) 

for functions u which belong to C 2 ((xi,X2))- Then for each u\ G \u(xi),u(xi)], U2 E 
[^(^1)) ^(^2)] there exists a unique solution u to (10. 1|) in [x\,X2\ satisfying u < u < 

U, u(xi) = Itl, u(x 2 ) = U 2 - 

This result is well known in several versions. It is proved, for example, in [S] in the 
context of ODEs and in [6] in the context of elliptic PDEs. 

Proposition 1.2. (i) Let xq € (0,oo). Suppose that there exist u,u G C°[0, xq] nC 2 (0, xq) 
satisfying 

n(0) = u(0) = (1.5) 

in addition to (l.l)-(l.S) with x\ = 0,^2 = xo Then, for each uq G \u(xo), u(xo)] there 
exists a solution of (0.1), (0.2) in [0,rco] such that u<u<uin (0, xq) and u(xq) = uq. 
(ii) Let u,u satisfy (l.l)-(l.S) for x\ = and X2 = 00 as well as (11 .5f) . Then, there exists 
a solution of (0.1), (0.2) in [0, 00) such that u< u<u 

Proof. (i) By Proposition 11.11 for each e G (0,xo) and each uq G \u(xo) ,u(xo)] 
there exists a unique solution u £ G C 2 ([e, xo]) of 

f £u £ = in (e,x ), 

\ u £ {e) = u(e), u £ {x ) = u , 

which satisfies 

u(x) < u £ (x) < u(x) for all x G (e, xq). (1-7) 

Combining these two-sided estimates with well-known elliptic regularity theory applied to 
(|1.6p . for each fixed 5 G (0, xo) we can find K(5) > such that 

\K\\cms,x }) < K{5) for all e G (0,<5). 
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As in the proof of Theorem 2.2 in [3], a straightforward extraction procedure involving 
the Ascoli-Arzela theorem enables us to pick a sequence (£k)keN C (0, Xo), u £k such that 
£k — > and u £k — > u in C^ c ((0, xo]) as /e — )• oo with some u 6 C 2 ((0,xo]) which, by (jl.6p . 
satisfies £u = in (0, xq). Furthermore, trivially u(xq) = uq and from fjl -T[) we obtain 
u < u < u in (0, rr ), hence 

= liminf < liminf u(x) < limsup u(a;) < lim sup u(:r) = 0, 
x-5>0 z->0 ^.^o 

(ii) By (i), for each e we have a solution of (0.1), (0.2) such that u £ (0) = and u<u<u 
in [0, 1/e]. Applying for e > the same extraction idea as in (i) we obtain the claimed 
solution in [0, oo). //// 



Proposition 1.3. For a + b > there is at least one solution of W.l\) , W.2\) on [0, oo) 
such that < u < x. 

Proof. For a + b = such a solution is u(x) = x. Suppose a + b > 0. Then, it 
can be readily checked that u(x) = is a subsolution and u{x) = x is a supersolution in 
[0, oo). The claim thus follows from Proposition 11.21 //// 



The problem (jO.ip . f|0.2|) can for a + b > be formally solved by a power series. We let 



and inductively define 



where 



for n > 1. 



fco:=l, ^ : =-^V^' (L8) 



/ n (s) := 5>z 1+i/2 , (1.9) 

i=0 

Cn+1 := 3ch(n + 3)x("+W (L10) 



Lemma 1.4. Xei a + 6 > 0. T/ien the coefficients {A;j}™ =0 are well-defined and £f n = 
0(x ( - n+2 ^ 2 ) as x \ /or a/2 n G N. 

Proof. The statement clearly holds for n = 1. Arguing by induction, we suppose 
that it is valid for some n > 1. Then 

£/n+i = £/„ " c (2(/; - 1)(1 + (n + l)/2)A: n+1 x("+ 1 )/ 2 + (1 + (n + l)/2) 2 k 2 n+1 x n+1 

= £f n -^ck 1 (n + 3)k n+1 x^ n+2 ^ 2 + 0(x^/ 2 ) asx\0. (1.11) 
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Since £f n is a polynomial in powers of yfx and £f n = 0(x^ n+2 ^ 2 ) it follows that k n +i is 
well defined and that Sf n - (3/2)cfci(n + ?,)k n+1 x^ n+2 ^ 2 = 0(x^ n+3 ^ 2 ). In view of (fTTTjl 
this implies that £f n +i = 0(x( n+3 )/ 2 ) and thus completes the proof. //// 



Note that the radius of convergence of the power series (jl.9p is nil and hence the power 
series does not define a solution directly via lim n _ 5 . 00 f n (x) outside x = 0. 

We will show later (Corollary 13. 4p that every solution u of (0.1), (0.2) satisfies 

n ( fc )( x ) = /W( x ) + ( x ( n+3 )/ 2 - fc ), 

for /c G {0, 1} and n = 1. Whether this is true for n > 1 or fc > 1 remains an open 
question. 

2 Nonexistence for a + 6 < 

In this second part we shall deduce Proposition 12.21 below which will exclude the existence 
of any continuous solution to (jO.ip for any xq > under the assumption a + b < which 
is complementary to the hypothesis of Proposition 11.31 

To this end we first prove that any supposedly existing continuous solution must satisfy 
u'(x) — > 1 as x — > 0. This property can formally easily be guessed upon tracing the 
possible solution behavior near x = 0. 

Lemma 2.1. Suppose that for some xq > 0, the function u E C ([0, xq\) n C 2 ((0, xq)) is 
a solution of W.l\) , W.S\) . Then 

lim tt'Or) = 1. (2.1) 

Proof. Letting v := u' — 1 we can rewrite (0.1) as 

u = v + 1, 
x 2 v = ax(v + 1) + bu — cv 2 . 

Let X(t) := -t' 1 for t < 0. Then X'(t) = t' 2 and X(t) \ as t -oo. We next 
introduce U(t) := u(X(t)) and V(t) := v(X(t)) for t < 0. Then the pair (U,V) solves the 
following system of differential equations 

u' = r 2 (v + i), 

V' = -af 1 ^ + 1) + bU - cV 2 . 
By assumption, we have U(t) — > as t — > — oo and thus 

V'(t)=p(t) + q(t)V(t)-cV 2 (t) (2.2) 

with 

p(t) -)• and g(t) -)• as t ->■ -oo. (2.3) 
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We wish to show that if V(t) is defined for all t < — Xq 1 then V(t) — > as t — > — oo. The 
proof proceeds in several steps. 

i) Given e > there is T < — Xq 1 such that \p(t)\ < e 2 c/3 and \q(t)\ < ec/3 for all t < T, 
by virtue of ([23]) . 

ii) Consider t < T. We claim that if | V(*)] > e for all i < t then 

nt) > yfa)-^ fa « S tc (2.4) 

while defined. To this end note that (j2.2[) and i) yield 

V'(t) < ~^V 2 (t) if t < T and |F(t)| > e. (2.5) 

By the comparison theorem for ordinary differential equations we conclude that 

V(t) > Y(t) for t < t , 

where Y solves the differential equation Y' = —^Y 2 with Y(to) = V(to). On solving for 
Y we obtain (f2T 



iii) Now we prove that there exists t± <T such that V(t±) > — e. Suppose to the contrary 
that V{t) < -e for all t < T. Then, (J23J) gives V{t) > -e/2 for t < t - ^, yielding the 
desired contradiction. 

Next we show that V(t) > — e for all t < t\. Arguing by contradiction, suppose this 
is not the case. Then there is t 2 such that — oo < t 2 = sup{t < t\ : V(t) < —e} < t\. By 
continuity we have Vfo) — ~ £ - From (|2.5p we obtain V'fo) < which is in conflict with 
V(t 2 ) = -e and V(t) > -e for t G (t 2 , h). 

iv) Finally, we show that V(t) < e for all t < T. If not, there is £3 < T such that 
^(£3) > e and we have t± := sup{i < : V(t) < e} < t 3 . The same argument as in part 
iii) shows that t± = —00 and therefore V(t) > e for all t <t%. From fj2.4j) we now obtain 
V(t) — > 00 for t \ £3 — |^(t3) -1 Therefore, V(i) is not defined for all t < — Xq 1 which is 
inconsistent with differentiability of U in (— 00, 0). 

Since e was arbitrary this completes the proof of the lemma. //// 



It is now possible to rule out local existence of a continuous solution of (jU.ip . (|U.2j) 
under the condition that a + b be strictly negative. 



Proposition 2.2. Suppose that a + b < 0. Then for each xq > 0, i/ie problem W.l\) , W.2\) 

does not possess any solution u € C ([0,xq]) fl C 2 ((0,xo)). 

Proof. Suppose that such a solution exists for some xo > 0. Then from Lemma 
12.11 we know that u actually belongs C 1 ([0,xo]) with u'(0) = 1, and hence the functions 
cpi and ip 2 defined by 

tpi(x) := v!{x) — 1, xG(0,xq), and ip 2 (x) := , xG(0,xq), 



7 



satisfy fi(x) — > and f2(x) — > as x — > 0. Since a + b < 0, we can thus find x G (0, xo) 
such that 

a + b + cupi (x) + 6^2 (^) < — ^— for all x G (0, x) . 

Therefore, (jO.ll) shows that 

x 2 u"(x) = axu'(x) + bu(x) — c(u'(x) — l) 2 

< axu'(x) + bu(x) 

= ax(l + <p\(x)) + bx(l + tp 2 (x)) 
a + b + a(pi (x) + b(f2 (x^j ■ x 

< —5x for all x G (0,x) 

holds with (5 := — > 0. By integration we find that 

~x 

u'Cx) — u(x) < —5 In— for all x G (0,x). 

x 

This implies that u'(x) — > +oo as x — > and thereby contradicts Lemma 12. 11 //// 



3 Monotonicity and concavity properties of solutions 

In this section we assume a + b > and we study monotonicity and convexity properties 
of solutions to (jO.ip . whose existence was established in Section 1. 

The following lemma is the key to establishing monotonicity, concavity, and ultimately 
also uniqueness in a certain restricted class of solutions. 

Lemma 3.1. Consider a nonconstant function y G C°([0, oo)) n C 2 ((0, oo)) satisfying 

x 2 y"(x) = f(x)y'(x) + g(x, y{x)), (3.1) 

for some continuous functions f and g. Suppose there is a constant y* G [—00,00] such 
that for all x > one has g(x,y) > for y > y* and g(x,y) < for y < y* . Then there 
is at most one xq G (0,oo) such that y'(xo) = 0. If such xq exists then one, and only one, 
of the following two alternatives is possible: Either 

y'(x) < for x < xq, y'(x) > for x > xq, and y(x) > y(xo) > y* for all x 7^ xq, 

or 

y'(x) > for x < xq, y'(x) < for x > xq, and y(x) < y(xo) < y* for all x 7^ xq. 
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Proof. We first note that because of continuity of g we have g(x,y*) = for all 
x > whenever y* is finite. By an ODE uniqueness argument, y(cco) = V* an d y'( x o) = 
implies y(x) = y* . Therefore, if y(x) is not constant and y'(xo) = then y(xo) ^ y* . 
Now suppose that y(xo) > y*. Then from (|3.1[) it follows that y"(xo) > 0, hence y'(x) < 
for x < xq sufficiently close to xq. Arguing by contradiction, let us suppose that there 
exists < x\ < xq such that y'{x\) > 0. Then there is X2 € [x±,Xq) such that 

y'{x 2 ) = 0, y'(x) < for x 2 < x < x , (3.2) 

which implies y"{x 2 ) > and also y(x) > y{xo) > y* for x 2 < x < Xq. On the other hand, 
(|3.ip together with y{x 2 ) > y* and y'(x 2 ) = entails that y"(x 2 ) < 0, yielding the desired 
contradiction. Therefore, y'(x) < for all x S (0,rco). The proof of y'{x) > for x > xq 
follows the same lines. Finally, the case y{x) < y* can be reduced to the case y(x) > y* 
by the transformation y \— > —y, y* i— > —y*. //// 

We now apply this to derive some monotonicity properties of solutions. Here in order 
to abbreviate notation, we call a function </> : [0, oo) — > K eventually monotonic if it is 
monotonic on [xo,oo) for some xq > 0. 

Lemma 3.2. Suppose that u G C 2 ((0,oo)) is a bounded nonconstant solution of IIP. 1]) on 
(0,oo). 

i) If b ^ and u is eventually monotonic, then u(x) converges to the unique stationary 
solution u = c/b as x — > oo. 7/6 = and u is eventually monotonic, then u is unbounded. 

ii) If b > 0, then one of the following alternatives occurs: Either 

u'(x) < for all x > 0, u(x) > c/b for all x > and u(x) — > c/b as x — > oo, 

or 

u'(x) > for all x > 0, and either u(x) < c/b for all x > and u{x) — > c/b as x — >■ oo, 
or u is unbounded, 

or finally 

there exists a unique xq > such that u'(xq) = 0, and we have 

u"(xo) > 0, u(x) > u(xo) > c/b for all x / xo, u'{x) < for x < xo, u'(x) > for x > xq, 
and u is unbounded. 

Hi) If b < then u'(x) > for all x > and u is unbounded. 

Proof. i) The substitution x(t) = e*, u(t) = u(x(t)) transforms equation (jO.ip into 

u" = (a + l)u + bu- c{e- l u - l) 2 . (3.3) 

Being bounded and eventually monotonic, u has a limit / as t — > oo and consequently 
lim^oo u'(t) = 0. From (j3.3H it now follows that lim^oo u"(t) = bl — c. If bl — c ^ then 
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lim^oo u"(t) 7^ which is inconsistent with the convergence of v! . This proves bl — c = 0. 
For 6 = this is a contradiction with c > 0, for b ^ it yields I = c/b. 

ii) If u'(x) < for all x then u is bounded and, by i), tends to c/b as x — > oo which is 
possible only if u(x) > c/b for all x. If u'{x) > for all x then thanks to monotonicity, 
u(x) approaches a limit in [0, oo] as x — > oo. If this limit is finite it has to equal c/b by 
virtue of i), and in the remaining case u is unbounded. 

Suppose now there is xq > such that u'{xq) = 0. Lemma I3TT1 applied to equation (|0.ip 
with y = u > 0, g(x, y) := by — c and y* := c/b yields two alternatives, the first of which is 
stated in part ii). The second alternative is not possible since it implies < u < c/b but 
at the same time u'(x) < for x > xq which means that u -ft c/b as x — > oo. A bounded 
solution not converging to c/b contradicts part i). 

iii) If v! is not positive everywhere then Lemma 13.11 applied to equation (jO.ip with 
y = u > 0, g(x, y) := by — c and y* := oo implies that there is xo such that u(x) < u(xo) 
and u'(x) < for x > xq. Therefore u is bounded and eventually monotonic. This contra- 
dicts i) when b = 0. For b < 0, i) dictates that u should converge to c/b as x — > oo, which 
contradicts u > since c/b < 0. //// 



Proposition 3.3. Suppose that a + b > and that u is a solution of (jO.ip . (|0.2p on 

(0,xq) with some xq > 0. Then there exists x\ S (0, xq) such that u'{x) > 0, u"{x) < 
and u"'(x) > for all x G (0, xi). Furthermore, in this case we have 



Proof. We recall that by Lemma 12.11 u'(0) = 1. As an immediate consequence we 
must have u'(x) > for all sufficiently small x > 0. On differentiating the equation (|0.ip 
we obtain 



Lemma I37L1 applied to equation (13. 5p with y = u' , g(x,y) := (a + b)y, y* := implies that 
u" has a constant non-zero sign near x = 0. This, together with u(x) < x and u'(0) = 1, 
yields that necessarily u"(x) < for all sufficiently small x > 0. 
We now differentiate equation (|3.5[) once more to obtain 

x 2 u"" = ({a - 4)x - 2c(v! - l))v!" + {2a + b- 2)u" - 2c{u"f on (0, x Q ). (3.6) 

Lemma 13.11 applied to equation (|3.6p with y = u" < 0, g(x, y) := (2a + b — 2)y — 2cy 2 
and y* := (2a + b — 2)/(2c) implies that u"'(x) has a constant non-zero sign near x = 0. 
Arguing by contradiction, we suppose that u'" < near x = 0. Since n" < 0, this 
implies that L := lim^o u"(x) exists and is finite. This however contradicts equation 
(10. ip . since on integrating we find x 2 u"(x) = Lx 2 + o(x 2 ), xu'(x) = x + Lx 2 + o(x 2 ) and 
u(x) = x + Lx 2 /2 + o(x 2 ) as x — > 0, and on substituting these expressions into equation 




(3.4) 



x 2 v!" + 2xu" = axu" + (a + 6)u' — 2c(u — l)u' 



on (0,x ). 



(3.5) 
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(|0.ip one concludes that it cannot hold near x = 0. We have thus proved u > near 
zero. 

Next, dividing (|0.ip by x we obtain 

xu"(x) = au'(x) + b + c ~ 1)2 for all x G (0,x ). (3.7) 

X X 

Since u" is increasing and negative, by (|2.ip we find that 

u'(x) - 1 = / u"(0d£ < < 0, 

and, consequently, 



o 



xu"(x) ^ as x -> 0. (3.8) 
Substituting this into f|3.7j> we obtain 

v (u'(x) - l) 2 

hm c v v ; — = a + b. 

x^o x 

Since u'(x) — 1 < 0, this is equivalent to (|3.4|) . //// 



Corollary 3.4. There is a continuum of local solutions of (0.1), (0.2), and they all satisfy 



u{x) = x--J ^x 3 / 2 + o(x 3 / 2 ) (3.9) 
3 V c 

u ' {x) = i _ ^/^x 1 / 2 + o^ 1 / 2 ). (3.10) 

PROOF. (jSHO]) follows from ([HD , and (H2D follows by integration of (JOU). //// 



Proposition 3.5. Let u be a solution of (jO.ip . (I0.2p wit/i xo = oo suc/j i/iai < n(x) < x 
/or a// x > 0. Then, in addition to (|3.9|) and (|3.10p . we /iave that u'{x) > 0, u"(x) < 
and u"'(x) > /or all x > 0. Moreover, 

lim n'(x) = 0. (3.11) 

Proof. The conclusion u'{x) > for all x > is a trivial consequence of Lemma 
13.21 iii) for b < 0. In the case b > 0, Lemma 13.21 ii) implies that if u is not increasing 
everywhere then there is xo > such that u'(x) < on (0, xq) and this contradicts the 
facts that u(0) = and u > 0. 

Next, Lemma [3J] applied to equation (13.51) with y = u' > 0, g(x,y) := {a + b)y > and 
y* := shows that if u" < does not hold over (0,oo) then u"(x) > for all sufficiently 
large x > 0. We show that the latter alternative is impossible. To this end, we let 
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v := v! and v(t) := v(x(t)) with x(t) = e t . Arguing by contradiction, since v is eventually 
increasing and u(x) < x, we must have t> < 1, which implies that v and 5 converge as 
x — > oo and i — )• oo, respectively. This in turn implies v'(x(t)) — > and v'(t) — > as 
t — > oo. Since from (|3.5p we see that 

5"(f) = ( a - 1)5' (i) + (a + 6)t;(t) - 2c(5(i) - l)v'(x(t)), (3.12) 

we therefore obtain v"(t) — (a + b)v(t) — > as i — >• oo. Since 5 > is increasing and 
convergent, so is v" . This is inconsistent with v'it) — > as t — > oo. We have thus proved 
< for all x > 0. 

Now since u is decreasing and bounded below by 0, it converges as t — > oo. Therefore, 
both v'(x(t)) and v'(t) converge to as t — > oo. From (|3.12p we thus obtain v"(t) — (a + 
b)v(t) — > as t — > oo. Should the limit of v for t — > oo not be zero, the same would hold 
for v" . This, however, is in conflict with the convergence of v' for t —> oo. This proves 
statement (13. lip . 

Finally, in order to verify that v!" > throughout (0, oo), we note that from Lemma 
13.31 we know that v!" is positive near zero. Arguing by contradiction, we assume that u" 
is not increasing everywhere. Lemma 13.11 applied to equation (j3.6[) implies that in such 
case u" must be eventually decreasing and therefore linx r _ >00 u"(x) ^ which contradicts 
statement (fBTTTT) . //// 



The following inequality is related to the theory of speculative attacks, in which an agency 
artificially supports a (low) fixed price of an asset, using a limited amount of reserves. The 
inequality indicates that in order for the speculator to make expected profit or at least for 
her to break even, the price of the asset must always jump upwards after the speculative 
attack has exhausted the entire supply of the asset at the subsidized price. 

(10. 2p with a + b> such that < u(x) < x 

for all x>0. (3.13) 

Proof. Since u"'(x) > for x > by Proposition 13.51 the function u' is strictly 
convex on [0, oo). Therefore 

u'(x) - u'(0) ,,, . 

— — — < u (x) for all x > 0. 

x 

On multiplying both sides by x, utilizing u'(0) = 1 and integrating we obtain 

u(x) — x < xu'(x) — u(x) for all x > 0, 
which yields the desired inequality. //// 



Corollary 3.6. Let u be a solution of (10. ip . 
for all x > 0. Then 

,, , u(x) 
l + u'(x) >2^- 

X 
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4 Uniqueness of the global solution bounded between 
and x 



Proposition 4.1. There is one, and only one, solution u of (jO.ip . (|0.2p on (0, oo) which 
has the additional property that < u(x) < x for all x > 0. This solution necessarily 
satisfies u > 0, u' > 0, u" < 0, and v!" > on (0, oo). 

The proposition stems from the following result: 

Lemma 4.2. Let u ^ v be two solutions of (|0.ip on (0, oo) which satisfy ii(0) = v(0), 
u'(0) = v'(0) and u" < on (0, oo). Then w := v — u satisfies either w" > on (0, oo) or 
w" < throughout (0,oo). 

Proof. The function w solves 

x 2 w" = axw' + bw — 2c(u' — l)w' — cw' 2 on (0, oo), 

which on differentiation yields 

x 2 w"' = {a- 2)xw" + (a + b)w' - leu 1 ' vJ - 2c(u' - l)w" - 2cw'w" on (0, oo). (4.1) 

Lemma [3.11 applied to equation (|4,ip with y = w' , g(x, y) := (a + b — 2cu" (x))y and y* := 
shows that w" can have at most one root. Now the existence of such a root xq > would 
imply either w' > 0, w"(x) < for < x < xq and w"{x) > for x > xq; or w' < 0, 
w"{x) > for < x < xq and w"{x) < for x > xq. This however contradicts w'(0) = 0. 
Thus one must have either w"(x) > or w"(x) < for all x > 0. //// 

Proof of Proposition 14.11 Suppose there are two solutions u ^ v bounded between 
and x and let w := v — u. Proposition 13.51 yields lim^oo u'{x) = lim^oo v'(x) = and 
therefore 

lim w'(x) = 0. (4.2) 

x— >oo 

On the other hand, Proposition 13.51 also gives u'(0) = v'(0) = 1, implying w/(0) = 0. We 
can thus employ Lemma 14.21 to obtain that either w"(x) > or w"(x) < for all x > 0. 
In view of w'(0) = both alternatives contradict (|4.2p . The claimed further properties of 
the unique solution bounded between and x follow from Proposition 13.51 //// 



5 Finer aspects of non-uniqueness 

By Corollary 13.41 there is a continuum of local solutions sharing the same asymptotics 
up to degree 3/2 at 0. The solutions are parametrized by their values at any x > 0. In 
this section we establish existence of disjoint continua of local solutions of (0.1), (0.2) 
distinguished by refinement of their asymtotics at 0. 
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We first show that the higher order terms specification of (|3.10p can be somewhat 
sharpened. Denote 

Lemma 5.1. Let u,xq,xi be as in Corollary \3.4\ Then there exists I > such that 

u'{x) - 1 



+ k 



<l\fx for all x & (0,xi). 



(5.1) 



If u is the unique solution provided by Proposition \4-l\ then (|5.ip holds for all < x < oo 
for some I > 0. 



Proof. Since 



such that 



Therefore trivially 



and also 



u'{x)-l 



-k as x \, by Proposition [331 we can find Xi > 0, c\ > 



bu{x) 


-b 






X 









|?/(x) — 1| < ci\/x for < x < xi. 
|a?/(x) — a\ < c%\a\y/x for all < x < x\ 

2ci|6| 



(5.2) 
(5.3) 



- l)dy 



< 



ci|6| 
x Jo 



\fydy = yfx for all < x < x\. 

(5.4) 

Taking x\ sufficiently small, using Proposition 13.31 we obtain that u" < and u'" > and 
hence \xu"(x)\ < \u'(x) — 1| for all < x < X\, from (15.21) we also infer that 

\xu"(x)\ < c\yfx for all < x < x\. 
Accordingly, from (jO.ip we obtain that (p(x) := - , < x < x\, satisfies 



cip 2 {x) — ck' 



(au'(x) -a)+ -bj- xu"(x) 



< C2y/x for all < x < x\ 



(5.5) 



with C2 := c\\a\ + Cl J 1 + c\. Therefore <p 2 {x) < k 2 + so that 



<^(x) = |^(x)| < A; • \l 1 + ^ < k ■ ( 1 + C2 ^ 



2dc 2 



for all < x < xi, (5.6) 



where we have used that vT+T < 1 + f for all £ > 0. Likewise, (|5.5j) entails that 
V? 2 (x) > k 2 - f for all s > 0. Thus, since y 7 ! 77 ! > 1 - ^ for all £ £ (0, |), we see 
that 

for all x E (0,x 2 ), (5.7) 



¥>(x) > A; - a/ 1-^5 > 



c/c 2 



4cfc 2 
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where X2 '■= min{xi, (§|j) 2 }. 

If u is the global solution of Proposition 14. 1\ then (|5.ip is obvious in [xi,oo) because 
u'(x) — > as x — > oo. This, together with (|5.6p and (J5.7P concludes the proof of the 
lemma. //// 



The core of our approach will be formed by the usage of on ordered pair of sub- and 
supersolutions which deviate from our original solution by an exponentially small term 
near x = 0. As a preparation, let us compute the action of the operator £, as defined in 
(|1.4p . on such functions. 

Lemma 5.2. Suppose that u is a solution of \0.1\) , in (0,xo) for some xq > 0. Then 

for all e > 0, a£l and E R, the function v defined by 



v{x) := u(x) — ex a e ^ , x £ (0, xq), 



satisfies 



£v 



ee ^ 



P 2 a -i 

+ 
+ 



a 



3\/3 + a/3 a/3 



2/ 2 



X 2 



a(a — 1) — aa — b 



■ x 



+2c- 



,a-l 



-ec ■ 



u'-l vP r 
—x 2 + ax 



+ ax 2 



2 _J_ 

■ e v^s 



(5.8) 



(5.9) 



for all x £ (0, xo). 

Proof. We write 

__§_ 

w[x) := ex a e , iG (0,so), 

so that v = u — w and consequently 

£ v = —x 2 u" + axu + on 
+x 2 w" — axw' — bw 
+c(u'-w'-l) 2 in(0,xo). 

Since (u' — w' — l) 2 = (u' — l) 2 — 2{u' — l)w' + w' 2 , using that £u = we see that 

£v = x 2 w" — axw' — bw + 2c(u — l)w' — cw' 2 in (0, xq). (5.10) 

We now compute 

w'(x)=e—x a 2e ^ + eax a e ^ 
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and 



w"(x) = £-^x a 3 e ^+e{a-^j^x a a e ^ + £^-x a a e + ea(a - l)x a 2 e 
and hence obtain from (|5.10p that 



Sv 



ee V s 



„a-l 



a/3 



+ ( a - -J |x Q "2 + -^-x a ~2 + a( a - l)x a 



P a — — a i a 

—x 2 — aax — ox 



+2c- 



—EC 



u'-l 



Pa- 



—x 2 + ax 



a-l 



Pa 



x 2 + ax 



a-l 



2 §_ 

e 



in (0,x ). 



On straightforward rearrangements, this yields (|5.9p . 



//// 



We can now identify an appropriate family of subsolutions for (jO.ip . 

Lemma 5.3. Suppose that for some xq > 0, u is a solution of \0.1\) , W.2\) in (0,xo) 
satisfying < f5. 1\) with some I > 0. Then for any real number a < | + a — 2cl there exists 
X3 = xs(a) G (0,xq) such that for each e > 0, the function v defined by 

4ck 

v(x) := u(x) — ex a e xG(0,x 3 ), (5.11) 

satisfies 

£v<0 in(0,x 3 ). (5.12) 
Proof. Given a < | + o — 2cl, thanks to the positivity of c and k we can find 



x\ G (0, xq) such that both 



2ck 



a > 



and 



2cA; • (| + a - 2cl - a) 



> a(a — 1) — aa — b + 2da 



(5.13) 



(5.14) 



hold. Then for e > we define v as in (15.111) and thus obtain from Lemma 15.21 applied 
to f3 := 4c/c, that (|5,9p holds for t; = v and all x G (0, x\). In order to utilize (|5.ip 
appropriately, we note that in view of (|5.13p we have 

-af-i + ax -* = (-= + a) ■ x a ~* >0 for all x G (0,^). 
2 V x/x / 



We may thus multiply this by the inequality 
u'(x) - 1 



<—k + l-/x for all x G (0, xi), 
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as resulting from (|5.ip . to infer from (|5.9j) on dropping a nonpositive term that 

?2 



Sv 



f> \ B 2 



ee ■ < —x 
I 4 

+ 
+ 



3\ (3 a/3 a/3 
° ~ 2/ 2 + ~2 2 



x a ~2 



a(a — 1) — aa — b 



_j_ 

ee V* 



+2c(-k + ly/x~) 
7? 2 



—x + ax 2 



L 4 

+ 

+ 



c/e/3 



a-1 



ol — — \ — h — — — 2c/ca + clB 
2/ 2 2 2 M . 



X 2 



a (a — 1) — aa — b + 2cZa 



in (0,xi) 



Here the leading term disappears because (3 = Ack, whereas 
3\8 a(3 a(3 „ , ,„ „ , r 



a I — I — : : 2cka + clB = 2ck ■ 

2/2 2 2 

= 2ck ■ 



a \- a — a — a + 2cl 

2 



a a + 2cl 

2 



< 0. 



Hence, using (|5.14p we infer that 



„ f 2ck ■ (4 + a - 2d - a) 
Sv < ee • x a ■ { ^ +a(a-l)-aa-b + 2cl 



< in (0,xi), 



as desired. 



Our construction of supersolutions to (10. If) is similar. 



//// 



Lemma 5.4. Let u be a solution of (0. W.Sfy in (0, xq) for some x$ > 0, which satisfies 
( 15. 1\) with some I > 0. Then for each eo > and a > | + a + 2cZ one can picfc £4 = 
£4(0) S (0, xq) such that for any e G (0, £o)> ^ e function v defined by 



satisfies 



v(x) := u(x) — ex a e , 3; G (0,14), 
£v > m (0, X4). 



(5.15) 
(5.16) 
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Proof. Relying on our assumption on a, let us first pick X4 G (0,xq) such that 
besides again 

^L+a>0, (5.17) 



V x 3 

the inequality 

2ck(a - § - a - 2cl) . . , . f „ r , „ 3 „ ,n2 _4c*f 

(5.18) 

holds. Then writing /3 := Ack, we again have fx a_1 + a/"2 > in (0, X4). Hence, using 



> |a(a— 1)— act— b— 2cla\+c- sup <J x c 



that (|5.ip implies that 

u'(z) - 1 



> —A; — /-v/x for all x G (0, Xq), 



from (15.91) we obtain 



frJ > ee 7? J ^Z" 1 



3\/3 a/3 a/3 
° ~ 27 2 + ~2 2~ 



x 2 



a(a — 1) — aa — b 



■ x 



+2c-(-k - ly/x) ■ 



-e c- 



P a-S , 

— x 2 + ax 



a-l 



e -S* 



in (0, X4) 



Since /3 = 4ck, this reduces to 

' 2ck(a - § - a - 2c/) 



4cfc 

£v > ee • x c 



\fx 



+ a(a — 1) — aa — 6 — 2cla 



-SqCX 



2ckx a 2 + ax 



a-l 



2 4cfe 



e v in (0, X4), 



and hence (|5.18p asserts ()5.16p . 



//// 



By means of Proposition 11.21 we can now infer the existence of infinitely many classes of 
continua of local solutions to (jO.ip , (|0.2|) . Here we shall strongly rely on the fact that the 
numbers a in Lemma 15.31 and Lemma 15.41 can be chosen in such a way that the functions 
v and v defined in (|5.11|) and (|5.15p are ordered appropriately. 
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Proposition 5.5. Suppose that for some xq > 0, u is a solution of h0.1\) . W.Sfy in (0, Xq) 
satisfying 115. 1\) with some I > 0. Then for all a G (— oo,|+a— 2d) anda £ (|+a+2c/,oo) 
there exists xq G (0,xq) such that given e G (0, 1) and any uq satisfying 

( 4cfc 4cfc \ 

u(x) — ex-e ^ , it(x) — ex a e ^ a£ x = So> 

(E2P, ( [O.i?]) possesses a solution u on (0,xq) fulfilling 

4cfc 4cfc 

u(x) — ex— e ^ < u(x) < u(x) — ex a e for all x G (0, x ) (5.19) 

as well as u(xo) = uq. 

Proof. We invoke Lemma 15.31 and Lemma 15.41 with 8q := 1 to obtain X3 G (0,xo) 
and X4 G (0, xq) such that for any e G (0, 1), the functions v and v defined by v(x) := 

_ 4ck _ 4cfc 

u(x) — ex—e ,x £ (0, X3), and v(x) := u(x) — ex a e ^ , x G (0, X4), have the properties 
Ev_ < on (0,^3) and £v > on (0,3:4). Writing i?o := min{l, X3, X4}, we moreover see 
using a < a that u < u throughout (0, xo). Therefore the claim results upon an application 
of Proposition fOT i) . //// 
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